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In this work, tiie spin dynamics of a single electron under parametric modulation of a lateral quantum dot's 
electrostatic potential in the presence of spin-orbit coupling is investigated. Numerical and theoretical calcu- 
lations demonstrate that, by squeezing and/or moving the electron's wave function, spin rotations with Rabi 
frequencies on the order of tens of megahertz can be achieved with experimentally accessible parameters in 
both parabolic and square lateral quantum dots. Applications of parametric excitations for determining spin- 
orbit coupling parameters and for increasing the spin polarization in the electronic ground are demonstrated. 

PACS numbers: 



I. INTRODUCTION 

Single spin operations are an important component for possible realizations of quantum dot-based quantum computerlQ In 
principle, single spin manipulations can be performed using electron spin resonance (ESR) techniques with the application of 
oscillating magnetic fields. A recent experiment utilizing ESR methods was used to generate spin rotations in a single electron 
lateral quantum dot at low static magnetic field strengths, where the required frequencies were on the order of hundreds of mega- 
hertz. However, extending ESR methods to microwave frequencies at low temperatures while minimizing the accompanying AC 
electric field in the semiconductor quantum dots is experimentally challenging. 

Instead of considering the AC electric fields as something to be avoided, AC electric fields can be used to control the electron 
spin in a quantum dot. Modulations of the electrostatic potential can readily couple to the electronic degree of freedom in a 
semiconductor quantum dot and can be performed on the nanosecond timescald^, which is roughly the timescale associated 
with the Zeeman energy in static fields on the order of a Tesla for GaAs quantum dots. However, in order to translate electronic 
control into spin control, there must exist an additional coupling between the spin and electronic degrees of freedom. 

The two main methods of coupling the spin and electronic degrees of freedom are either through the Zeeman interaction or 

through spin-orbit coupling. Consider first the Zeeman interaction, which is given by Hz = B{r) ■ (^{rjfXBS^ , where B{r) is 

the applied magnetic field, is the effective g-tensor, is the Bohr magneton, and S is the electron's spin vector. Previous 
experimental work has demonstrated that electric fields can be used to move the electron around in the presence of a spatially 
dependent g-tensor, g(r), in order to generate an effective time-dependent Zeeman interaction, Hz{t) = h,Lij{t) ■ S, which was 
used to perform spin rotations^ ^. Alternatively, Hz{t) can also be generated by using an electric field to move an electron 
around in the presence of a spatially inhomogeneous magnetic field^. 

Besides manipulating the effective Zeeman interaction, modulating the electrostatic potential can be used to manipulate the 
electron's spin through the spin-orbit interaction. Previous theoretical studies have proposed using electric fields to generate spin 
rotations via an electric dipole spin resonance (EDSR); theoretical calculations indicate that EDSR methods can be used to rotate 
an electron's spin on the order of tens of nanoseconds for electrons in quantum well d'^Pi'" 'and quantum dots" For parabolic 
quantum dots, application of an electric field is equivalent to modulating the center of the quantum dot which, in the presence of 
spin-orbit coupling, is able to induce spin transitions when the modulation frequency is equal to the Zeeman frequency. Besides 
utilizing EDSR for spin excitation, electric fields can also be used to modulate the Rashba spin-orbit coupling constani^^in order 
to perform combined spin and orbital excitations in quantum dotP. 

In this work, the spin dynamics under the combined action of spin-orbit coupling and parametric modulations of the electro- 
static confining potential (such as squeezing and moving the electronic wave function) in both parabolic and square quantum 
dots is examined. Theoretical and numerical calculations utilizing Floquet theory and effective Hamiltonian theory are used to 
demonstrate that parametric modulations of the electrostatic potential can generate single spin rotations on the timescale of tens 
of nanoseconds using experimentally accessible parameters. In order to maximize the amplitude of the Rabi oscillations, the 
modulation frequency must be chosen with precision on the order of the Rabi frequency (order of 1-10 MHz). Such parametric 
excitations therefore can provide better energy resolution of the quantum dot's energy levels over transport measurements, where 
the measured energy levels are often poorly resolved due to th ermal effects. Besides performing single spin rotations, we have 
shown that measurement of the observed Rabi oscillationJISESIand the required modulation frequencies can be used to determine 
the spin-orbit coupling parameters in the quantum dot. Finally, we have demonstrated that squeezing and expanding the elec- 
tronic wave function in a parabolic quantum dot can induce both orbital and spin excitations, which can be used to increase the 
spin polarization of the lowest electronic state. It should be noted that relaxation is neglected throughout most of the paper since 
the calculated Rabi frequencies (on the order of 1-10 MHz) are at least one order of magnitude greater than the experimentally 
observed 1/Ti values found in GaAs quantum dots'^ i9 20^ 
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where 



II. PARAMETRIC MODULATIONS AND FLOQUET THEORY 

Before studying the case of parametric excitations in lateral quantum dots, the general formalism of using parametric mod- 
ulations to generate excitations in quantum systems is presented. Consider a Hamiltonian which is a function of the parameter 

uj: 

Hiuo) = V E . i-^) '^^^^^ ^^^^^ ' I V ^""^ '^'^^^^ ^^^^^ ' ^ (^) O'(^) I 

= W{uj)Ha{oj)W\uj) (1) 

5- ^ ^ p / J^o)(fco| . Ajfe(cj)|jo)(fco| + Afcj(w)|fco)Oo| 

^ ^ ^ |fcM)(fcol (3) 

with |fco) = |fc(cL'o)). The transformation W{uj) simply switches from the basis {|fc(ajo))} to the basis {|fc(a)))}. 

In order to induce transitions between two arbitrary states \k{iS)) and \i{Lo)) efficiently, the matrix element coupling these 
states must be modulated at a frequency ujr given by nhujr ~ \Ej{Ld) — Ek{Lo) \ for integer n. The time-dependent coupling 
between states can in principle be generated by modulating the parameter lo of H{oj). The propagator during modulation of lo 
can be written as: 

U{t) ^ Texp^-^^ dt'H{ijj{t')) 

= W{iu{t))Texp(-'- I dt'H{t')\ (4) 



where 



H(t) ^ H,{u:(t))-inw\u{t)f-^^ 



(5) 



with 



for j ^ k 



ih du;{t) (fc(a;ft))|g|gg»|j(^(t))) 



(7) 



for j ^ k 



where Ajfc(w(i)) — Ej{Lo{t)) — Ek{uj{t)). Note that J^,* dt'Vfefe(i')/?i is just the geometric phas^^Xl^ssociated with the parametric 
modulation of H{ui). 
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In the following, we are interested in the case where only small, periodic modulations to lo are performed, i.e, Lo{t) = 
L0 + 5bj sm{u)rt + 0) where 5u) ^ oj, and ujr and are the modulation frequency and the initial phase respectively. If the energy 
difference between the states |jo) and jfep) is much greater than the generated off-diagonal matrix elements (|Ajj;(a;(0))| ^ 
|Vj7c|), then the two states will remain uncoupled unless |Ajfc([x)(0))| is equal to a multiple of fiujr, i.e., nTiujr ~ l^jfc('^(0)li 
where n is an integer. In order to gain insight into the dynamics under periodic modulation of the Hamiltonian, Floquet theory 21 
can be used to solve for U{t) in Eq. (jijl. Rewriting H{t) in Eq. (jsjl as H{t) = Hq{ijj{Q)) + V{t), with i/o(w(0)) given by Eq. 
([sj and V{t) — Vm exp(imwri), the effective propagator can be written as: 

dt'H{t')^ (9) 
Hq{uj{Q)) + ^ Vm exp(imwrt) J 

Tn J / 

eyi]i(inuJrt)Un{t) (10) 




where Eq. ( 10 1 represents the Fourier decomposition of the propagator U (t) into the Fourier operators, Un{t), which satisfy: 



dt n 

m 

with the initial conditions chosen to be [/o(0) = 1 and C/„^o = in order that J7(0) = 1. The solution to Eq. 
|TT| i is given by U(t) = exp(— iiJi?t/?i)U(0), where U(t) can be formally represented as a vector in Floquet space, 
U(t) = Yl°^=-oc^n{t)\np) with U(0) = 1|0f) [where the various Floquet states \np) satisfy {np\mF) = 5nF,inF], and 
Un{t) — [np] exp (^—iHpt/hj \0p). The time-independent Floquet Hamiltonian, Hp, is defined by: 

{np\Hp\mp) = (np\HQ,p + NpTiLOrlmp) + (np\Vp\mp) 



(^Ho{uj{0)) + nphiOrl^ + J2 SuF-k^m^Vk (12) 



where the states in Floquet space are denoted by \j, np), with {Hq^p + NpTitUr) \j, np) = {Ej{uj{0)) + nFhUr)\j,np). 

By going into Floquet space, the initial time-dependent Hamiltonian, H{t), has been replaced by the time-independent Floquet 
Hamiltonian, Hp; however, in order to calculate U(t), we must now exponentiate Hp, which is an infinite dimensional matrix 
in Floquet space, since {np\Hp\mp) is defined for all n^-, G {— oo, oo}. However, the propagator can be simplified if there 
exists a finite dimensional subspace of nearly degenerate states, Q, which is weakly coupled to states outside of Q. In this case, 
the dynamics within Q can be separated from the rest of the Floquet space by constructing an effective Hamiltonian^^ within 
the subspace Q by treating the coupling to states outside of Q perturbatively. For example, take Q = {\j, np), \k, {n + p)p)} 
with Ej{uj{0)) ~ Ek{uj{0)) + pfiujr, and assume that Q is weakly coupled to states outside of Q (i.e., \ {j,np\Vp\q,mp)\ ^ 
\Ej{Lu{0))- Eg{uj{0)) + {n-Tn) phujrl and \{k, {n+p)p\Vp\q,mp) \ < \Et,{u}{0))-Eg{uj{0)) + {n+p-m)phuJr\), an effective 
Hamiltonian can be written in the subspace Q by using a unitary transformation, exp (^Sp^ , where exp (^Sp^ Hp exp (^—Sp^ = 

H^; the state \j, np) is only coupled to the state |fc, (n + p) p) in H^^, i.e., {p, rp\H^^\m, qp) ^ only if \m, qp) , \p, rp) e 
Q. A perturbation series for Sp in powers of Vp, Sp — X^i^i s'j^^ can be constructed and is given in Appendix A. 
Using H^^, the propagator, U{t) can finally be written as: 

oo 

W''{t)U{t) = exp{iLUrt)Un{t) 

n— OO 

oo / ■ \ 

— ^ exp(ina;rt)(n_p| exp (^—Sp^ exp (—^H^^tj exp (^Sp^ \0f) 

n— — oo ^ ^ 

« £ exp(ina;,0 ({np\ exp f-^Hf^t\ \0p) - {np0^\exp f-^Hf^XlOp)] 

n— — oo \ \ / \ / / 

(13) 
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The resulting propagator, projected onto the {|jo), l^o)} subspace (where = bo)Oo| + |^o)(fco|). is given by: 



where 



{j,np\VF\j.,OF) {j,(7l-p)p\VF\k,0) 

Y 1 nuJr Ajk-\-{n—p)huJr 

" ^ ' (fc,(«+p)f |V>|j,0) (fc,nF|V-f Ifc.Of) 

{n-\-p)huJr — Ajk nujj. 

Z{LUrt) = [I r (15) 



and U^^{t) — exp {^—iH^^t/fij, where Heff is th^ effective Hamiltonian in the {IjjOf), \k,PF)} Floquet subspace, i.e.. 

Hl^^ = {\],QF){3,QF\ + \k,PF){k,PF\)HfH\3,QF){3,QF\ + \k,PF){k,pF\). 



III. PARAMETRIC MODULATIONS IN LATERAL, PARABOLIC QUANTUM DOTS 

Using the formalism presented in Section |ll] we are now ready to begin studying parametric excitation in a single-electron, 
lateral quantum dot, which is taken to lie in the XY plane with the electron's wave function strongly confined along the z 
direction. In the following, only an in-plane magnetic field, B — Bcos{9)x + Bsm{d)y, will be considered, which allows one 
to neglect orbital effects associated with an out of plane magnetic field. The Hamiltonian for a single-electron lateral quantum 
dot defined by the electrostatic potential V{X, Y), in the presence of spin-orbit coupling (for the moment, only Rashbd^and 
linear Dresselhaus^S-are considered), is given by: 

(16) 

crz (17) 

{C2{-0)^x + Ci{0)^z) (18) 



H = 


Ho + Hso 




Ho = 


Px H 

+ + 


v{xS)- 


2m* 2m* 


Hso = 




+ (:2{e)dz) 



2 

Py 

h 



where Tiujz = |.9/^b^|> Ci(^') = acos{9) — (3sm{9), C,2{6) — asm{6) — (3cos{9) [where a and (3 are the Rashba and linear 
Dresselhaus coupling constants respectively], and the vector potential, A — Zsm{9)x — Z cos{6)y, was chosen. Due to the 



strong confinement along the z-direction, all terms hnear in Z have been truncated/removed from H in Eqs. ( 16 1-( 18 i, with terms 
quadratic in Z being incorporated into the the confining potential along the z-direction. The electrostatic potential, V{X,Y), 
mostly results from voltages applied to surface gates above the 2DEG, which confines the electron within the quantum dot; 
changing the voltages of the surface gates can change V{X,Y). The spin quantization axis has been taken to be along the 
direction of the in-plane magnetic field. The eigenstates of Hq are denoted by \n, ±), where Ho\n, ±) = {En T ^wz/2)|ri, ±). 
In the presence of spin-orbit coupling, the various \n, ±) states are mixed; however, if the confinement strength is much larger 
than the spin-orbit coupling strength, i.e., \ {n,±\Hso\Tn,±)\ <C |A„„i| and \ {n,^\Hso\^,^) \ |A„„i±wz| for all n and 
m, Hso is suppressed and can be treated as a perturbation to Hq- 

In this section, a lateral quantum dot defined by a parabolic electrostatic potential will be examined. Such parabolic potentials 
have enjoyed tremendous success in describing transport and spectral properties in lateral quantum dots^^ The electrostatic 
potential for a parabolic quantum dot can be written as: 

nX,Y) - ^{x~.,)\!I^(Y-y,y^eF.X-eFyY 

m*ujj^ J 2 \ m*ujY ^ / 2m*ujj^ 2m*ujY 

(19) 



where = Xc'X + ycS is the center of the parabolic well, u)x and ijjy are the oscillator frequencies of the parabolic 
dot, and Fx and Fy are static electric fields which do not alter the energy levels of the quantum dot but do shift the ef- 
fective center of the quantum dot to Tc = [xd + eFx/(m*a;^)]x + [yc' + eFy/(m*ajy)]y. For convenience, V{X,Y) 
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FIG. 1: Parametric modulations of a parabolic quantum dot's confining potential by either modulating (A) the confining frequency (ljx and/or 
ljy) or (B) the center of the dot, fc — XcX + i/cy. In the presence of spin-orbit coupling, if either Sujx(y) (*) or Sfc{t) = 5xc{t)x + Syc(t)y 
is modulated at roughly the energy difference between |0, 0, +) and |0, 0, — ), effective spin rotations can be performed as depicted in Fig.[T| 



in Eq. (19i was taken to be separable in the X and Y degrees of freedom; that is, the principal axes of V{X,Y) were 



chosen for convenience to be the same as those used in Eq. (16i (if this is not the case, inclusion of terms like \XY 
can be readily incorporated into the following theory). The eigenstates of for parabolic confinement,|ri, m, ±), satisfy 
HQ\n^m,±:) = [hLOx{n + \/2) + huiYim + \ /2) ^fibJz /'i\\n,ra^±) and are centered about Tc. Under conditions of strong 
confinement (where TiijJx , Ttfj^Y ^ | i?5o I )> the states |0,0,+) and |0,0,— ) are approximately the two lowest energy eigenstates. 
Transitions between these two states correspond to spin rotations within the ground electronic state. 

As discussed in Section|ll] efficient transitions between the states |0, 0, +) and |0, 0, — ) in the quantum dot can occur whenever 
the Hamiltonian is parametrically modulated at roughly the energy difference between the two states. There are two natural 
modulation parameters to choose from in parabolic quantum dot potentials as shown in Fig. [T] either modulation of the oscillator 
frequencies {ojx(y) <^x(f) + ScuxiY) in Fig- [ijA)], or modulation of the center of the quantum dot [rc rc + 5fc in Fig. 
[TJB)] (note that modulation of Fx and/or Fy is equivalent to modulating r^,). Such parametric modulations can in principle 
be performed by modulating the surface gate voltages which define V{X, Y). Modulations of surface gate voltages have been 
experimentally performed on the nanosecond timescald^^. 



A. Parametric Modulations of the Confining Strength in a ParaboUc Quantum Dot 



The modulation the the oscillator strengths of the quantum dot, ujx(y) ^x{Y){t) = ^x(y) + Sujx(Y){t), corresponds to 
squeezing and expanding the electron's wave function in a time-dependent manner, as illustrated in Fig. [TjA). Using Eq. Q, the 
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effective time-dependent Hamiltonian during modulation of ojx and ojy at a frequency Ur is given by: 

^ = 0Jx{t){a^xax + ll'i) + uJY{t){a\.aY + l/2)-^dz 



' m 



.atJx(i) eFx / 2 / I \ .au;Y(i) eFy 



(20) 



= ^ ^go(tJx, ^y) + K» exp{imujrt)^ (21) 

where a^^y^ and ax(Y) are the creation and annihilation operators associated with the harmonic potential, and ujx(Y){t) = 

LUx(Y) + <^t^x(y) sin(cL'.rt + 0x(y)) is the time-dependent oscillator frequency. Note that Ho{u!x ,^y) contains the time- 
independent contributions to the energy from the harmonic potential and the Zeeman and spin-orbit interactions. In the presence 
of a static electric field (i.e., Fx{y)^ 0), modulations of uix and toy also results in a modulation of rc, which leads to the terms 



linear in ax(Y) and o^xiY) Eq. (20 1. 

rpToquet states |1) = |0,0,+,1f) and |2) = |0,0,— ,0f) can be performed if the modulation 
frequency, LOr, equals the energy difference between states |1) and |2), i.e., lo^ ~ ujz- Using Eq. Q, the effective Hamiltonian, 

^ ,2 

'^X{Y) 



C{Y 

Efficient transitions between the 
equency, ujr, equals the energy di 
i?EFF' ill presence of parametric modulation of the oscillator strengths at ojr ~ ujz is given [up to Sajj^^y)]- 



2 

where ^ ^^^^^^ _ \2){2\,d\^ = \l){2\,dl^ = |2)(l|,and 



-^EFF = ^{iOr-uJz- - 6z) d'i + b+d^^ + (22) 



^12 

by. - 



Tn*uz (LUx{Ci{^0)f , c^y(C2M))' 



m*ujz5ujl {C2{~0)f /32w8 + 4wf.a;| - 85wf,w| + 27wf-w| + 4w| 



+ 



^ V ^xi^z-^x) t^y(^l-^y) / 

^ eFxefy^LJx'^^y . , ^ v, /ly ^K^i + ^y) + ^i^y - 3^1 
2 -4 COS(0x - 0yjwzCi(~t'jC2(c'j 2^2/ 2 2T2 

f Ci{~(^)S^xeFx exp{±i(l)x) C2(-fi')<5wye-Fy exp(±i0y) ' 



V ^x (^i - w|) (wf, - a;|) 

.m*c^z / , ,Ci(-0)C2(^)'5c^x , exp(±i0y)CiWC2(-^)5^y 
± i_ (^exp(±i0x) + ^^^^^ 



(23) 



where A^^ is the higher-order contribution of the spin-orbit coupling to the energy difference (in Eq. (23 1, A^^ is only written 
to second-order in the spin-orbit coupling), and Sz is a Bloch-Siegert shift-^" which arises from the noncommutivity of H{t) at 



different times. In the presence of a static electric field, the coupling between the states |1) and |2), 6± in Eq. (23 1, is first-order 
in the spin-orbit coupling parameters, whereas it is second-order in the spin-orbit coupling when Fx — Fy — eV/m. 

Figure |2| A) presents the exact numerical simulation of the transition amplitude, |(0, 0, +|C/(t, 0)|0, 0, — ) p, for a parabolic 



quantum dot under harmonic modulation of ujx for Fx 7^ 0. In the simulation, the eigenstates of Ho{uJx ,ujy) in Eq. (20 1 were 
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numerically found by diagonalizing Hq{ujx ,^y) using a basis of four hundred |n, m, ±) states. Next, H{t) in Eq. (21 1 was 
expanded in the n lowest energy eigenstates of _ffo(^x, ^y), and U{t) = T exp{—i/1i J^^ H{t')dt') was then found numerically 
(n = 30 was found to give converged results for the simulations). The modulation frequency, uj^, was given by the numerically 
calculated energy difference between the two states, AE12, plus the Bloch- Siegert shift given in Eq. (23 1, i.e., cOr = AEi2/'h + 
5z ~ ijJz + A^z + ^z- This method was used throughout the paper when numerically calculating the exact transition amplitudes. 

The following parameters were used in the simulation shown in Fig. |2jA): Tlluy — T^'-^x — 1 meV, Tiujz = 0.1 meV 
(corresponding to an in-plane magnetic field of around five Tesla in GaAs), a = 4 x lO^^'^ eV-m, m* = 0.067too (where mo is 
the free electron mass), 6ujx — ^x/^^. Fx = lO** eV/m, and 6' = 0. With these parameters, the effective Rabi frequency for 
"on-resonance" modulation [ (w^ — W2)/(27r) = —3.46 MHz] observed in Fig. |2];A) was 9.86 MHz, which is consistent with 
the theoretical value given by Eq. ( [23] l, |(5± |/(27r) = 9.76 MHz. The fact that the Rabi frequency is directly proportional to the 
spin-orbit coupling arises because modulation of ojx results in a modulation of the effective center of the quantum dot Tc by dvc 
when Fx 7^ 0. For the parameters used in Fig. |2jA), \5tc\ = 2.3 nm. 

For the case when Fx — Fy — eV/m, remains fixed, and the electron wave function is only squeezed/expanded 
during the modulation [Fig. [ijA)]. Figure |2|B) presents the exact numerical simulation of the transition amplitude, 
|(0,0,+|L/(i,0)|0,0,-)p, when Fx = Fy = eV/m. The following parameters were used: huy = 1 meV, fiujx = 0.25 me V, 
fiLOz = 0.1 meV, a = 8 x lO"^'^ eV-m, m* = 0.067mo (where mo is the free electron mass), 5uox — t^x /20, and 6 — 7r/4. With 
these parameters, the effective Rabi frequency observed in Fig. |2];B) was 355 kHz, which is consistent with the calculated value 
given by Eq. (23 1, |(5±|/(27r) = 354 kHz, for "on-resonance" irradiation, (tjj, — (jj2)/(27r) — —39.285 MHz. The observed Rabi 
frequency was over an order of magnitude smaller than the case of nonzero Fx [Fig. |2|A)]. As stated earlier, the effective Rabi 
frequency is smaller since it is second-order in the spin-orbit coupling and contains terms like (a^ — ) /2 sin(26') ± a/3 cos(26') . 
Note that for 6 — mT/2, both a and j3 must be nonzero for 5± to be nonzero. Furthermore, for the case when ojx = i^Y and for 
uniform modulation (i.e, Sluy = 5uox), the Rabi frequency in Eq. (23 1 is exactly proportional to a(3 cos{26). Finally, it should 
be noted that oj^ — toz is on the order of |(5± | [Fig. |2|A)] or much greater than \6± \ [Fig. |2jB)], so that w,. must be tuned with 
precision given byJ5± | in order to maximize the amplitude of the Rabi oscillations. 

Although Eq. (mi was used to generate iJ^pp in Eqs. (22i-(23i, the dynamics under modulations of uix and can be 
calculated using the following time-dependent Hamiltonian: 



H{t) 



p2 
^X 



p2 



2m* 2m* 



Hso 



^z+'^Ux{t)X^ + ^Ut)Y' 



eFxX - eFyY 



(24) 



where Px{y) = y/m*huJx{Y)/'2Kax{Y) ~ "-xiY)), X{Y) = y^n/{2jn*ujx(Y)) (ax(Y) + a^jy^), and wx(y)(i) = ujx{y) + 
Sujx(Y){t)- H{t) in Eq. \2A\ can be used to construct an effective Hamiltonian in the |0, 0, +, If) and |0, 0, — , Of) subspace. 



which is also given by Eqs7122|-( 23 1 



B. Parametric Modulations of a Parabolic Quantum Dot's center 



In addition to modulating the confinement frequency, the center of the parabolic quantum dot can also be modulated, i.e, 
Xc + Sxc{t) and Uc ^ yc + Sydt). Unlike the case of modulating ujx and ojy, displacements of the harmonic potential do 
not alter the energy spacings of the quantum dot; however, the simultaneous eigenstates of the oscillator do change under such 
modulation. Using Eq. (j4|, the Hamiltonian under time-dependent displacements of the parabolic potential's center is given by: 

F[{f) 4. 4. ^Z--^ 
— ^ = uJx{a^x°'X + '^/'^) + ^^Y [ay ay + 1/2) —az 



■ / t \ / , / , /^N- N . Im*ujx dxJt) I 

(4 - «x) (Ci(-^)-x + c,moz) - ( 

• \'^*^Y / t \ (/■ (a\- , /■ f a\- \ ■ [rn*u}Ydydt) / 



a\ - ax 



ay (25) 



For small displacements, Xc ^ Xc + SxcSmioJrt + <t>x) and yc ^ Vc + SycSm{ujrt + (py), the effective Hamiltonian in the 
|0, 0, +, 1 f) and |0, 0, — , Of) subspace can be written as (for w,- ~ ujz)' 

rrl2 1 

= - - + 0^ + j 0-2 + o_|_o-_,_ + d_(T_ (26) 
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time (f^s) 

FIG. 2: Numerical simulation of the transition amplitudes, | (0, 0, +|rexp(— i dt' H{t'))\0,0, —)\^ , caused by modulating a quantum dot's 
oscillator frequency, ujx{t) ~ i^x + Sujx sin(ti;,.t), (A) with and (B) without a static electric field. In the presence of a static electric field, 
the effective Rabi frequency is first-order in the spin-orbit coupling [Eq. l|23[)], which lead to a large Rabi frequency as shown in Fig. |2jA). 
The following parameters were used in the simulation: Tiluy = ficux = 1 meV, a — 4 x 10~^^ eV-m, (3 = eV-m, huiz ~ 0.1 meV, 
9 = 0, Slox = Fx = eV/m, and (j>x = 0, which gave a Rabi frequency of 9.86 MHz compared with the calculated value of 
\5±\/{2n) ~ 9.77 MHz by Eq. \23\ . In the absence of a static electric field, the effective Rabi frequency is second-order in the spin-orbit 
coupling, which lead to a smaller Rabi frequency as shown in Fig. |2jB). The following parameters were used: hu>Y = 1 meV, hux = 0.25 
meV, hujz = 0.1 meV, a = 8 x 10"" eV-m, SuJx = lox /20, (3 = eV-m, 9 = 7r/4, Fx = eV/m, and (j^x = 0, which gave a Rabi 
frequency of 355 kHz compared with the calculated value of |(5±|/(27r) = 354 kHz given by Eq. \23\ . 



where 



+ ^ -J^COs(0S,-#)%fe,^zCl(-^)C2(e)^H 2 2r 2 2^2 (27) 

SI = [exp{±i<j)'x)Sxc 2 2 exp(±i'/>y)% 2 2 (28) 
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FIG. 3: Numerical simulation of the transition amplitudes, |{0, 0, + |rexp(— i dt' H(t'))\0,0, —)\^ , caused by modulating the center of the 
quantum dot's parabolic potential. In this case, the effective Rabi frequency is first-order in the spin-orbit coupling [Eq. ([28j], leading to a 
large Rabi frequency. The following parameters were used in the simulation: hcoY = hcox = 1 meV, a = 4 x 10~^^ eV-m, /3 — eV-m, 
TiUJz = 0.1 meV, 9 = 0, 5xc = 11.5 nm [equivalent to using an electric field of Ex = lO"* eV/m], and 4>x = 0, which gave a Rabi frequency 



of 48.7 MHz compared with calculated value of \S±\/{2n) ^ 48.8 MHz given by Eq. 1 28 



The effective coupling between spin states, S^., is first-order in the spin-orbit coupling, which leads to large Rabi frequencies. 

Figure Isjpresents an exact numerical simulation of the transition amplitude, | (0, 0, +|T exp(— i J* dt'H{t') /K)\0, 0, — ) p, in the 
presence of modulating the center of the parabolic well. The following parameters were used in the simulation: fiojx = Tt-loy = 1 
meV, a = 4 X IQ-^^ gV-m, /3 = eV-m, Uajz = 0.1 meV, 6* = 0, Fx = -Ff = eV/m, and 5xc = 11.4 nm. Such 
parameters gave an effective Rabi frequency of 48.7 MHz in Fig. Islcompared with the calculated value of |(5^j_|/(27r) = 48.8 
MHz given by Eq. ( 28 1 for "on-resonance" modulation (i.e., (tj^ —t^zj /(27r) = —3.435 MHz). As mentioned in the introduction, 
parametric modulation of a parabolic quantum dot's center is equivalent to applying a time-dependent electric field, E{t) — 
i?x(i)^ + -£'y(t)y, with the correspondence that 5a;c(t) = eEx{t)/{rn*Lo\) and Sydt) — eEyit) / [m*U)Y)- The simulation in 
Fig. Iljcorresponds to Ex = 10^ eV/m. Previous theoretical work has suggested using such EDSR methods to perform efficient 
spin rotations in quantum well&^-^-^and quantum dot a" '^ ' . 



C. Determining tlie spin-orbit coupling constants, a and /3 

In the above examples of parametric oscillations, both the effective Rabi frequency, and the effective offset, A^, depend 
upon the spin-orbit coupling constants of the quantum dot (a and (3), the oscillator frequencies {lox and uoy), and the direction 
of the magnetic field, 6. This dependence can potentially be used to determine both the spin-orbit coupling constants a and j3. 
For instance, consider the effects of modulating the center of the parabolic dot, i.e., performing an EDSR experiment. Taking 



the ratio of the "measured" Rabi frequency [Eq. (28 1] for an experiment with 5xc — XmoA and 5yc = to the "measured" Rabi 
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frequency for an experiment with 5xc = and 5yc = Amod gives 

d±{Sxc = Xmod, Syc = 0), 



ZZ{9) 



S±{6xc = 0, Syc = Amod, 



U-0) '''' 

For an in-plane magnetic field along the y axis (i.e., — 7r/2), the relative ratio of the spin-orbit coupling strengths can be found 
since in this case ZZ{tt/2) = [3 /a. Furthermore, A^, which can be determined by tuning LOr to the frequency which maximizes 
the amplitude of the Rabi oscillations and then approximating k, ujr — loz, can then be used to determine the absolute value 
of a and (3. Although the calculations used in this section were for a harmonic potential, more realistic V{X, Y) using the full 
electrostatic potential generated by the metallic surface gate could be used to more accurately characterize the spin dynamics 
in terms of a and /? under parametric modulation ofV{X,Y). 



IV. THE EFFECTS OF THE CUBIC DRESSELHAUS SPIN-ORBIT INTERACTION 

Recent work has indicated that the cubic Dresselhaus term can become the dominant spin-orbit interaction in highly confined 
(i.e., large lox and luy) quantum dots^^. For the coordinate system chosen in this work, the cubic Dresselhaus term is given by: 

H'd - ^(PYPxPY{com^z~^i^{0)dx)-PxPYPx{^HO)dz + cos{e)dx)) (30) 



The linear Dresselhaus coupling, /3 in Eq. ( [T8] l, is related to the cubic Dresselhaus coupling constant, 7 in Eq. (30i, by 
(3 ~ "f{Pz)/h^ [where for a quantum well of width (Pz) /^'^ ~ (tt/Zu,)^]. For a quantum well with 1^ ~ 30 nm, and 
for 'huJx(Y) ~ 1 meV, the relative strength of the cubic and linear Dresselhaus terms is given roughly by {Px(y)) I O^z) ~ 
TO*wx(y)^?u/(27r'^7i.) — 0.04. The larger ujx and uuy are for the quantum dot, the more the cubic Dresselhaus coupling con- 
tributes to the total spin-orbit interaction. 

Under parametric modulation of the center of the quantum dot (i.e., an EDSR experiment), the effect of the cubic Dresselhaus 
coupling can be found by simply adding i?™'' to H{i) in Eq. (25i. The effective Hamiltonian in the |1) = |0,0,+, \f) and 
|2) = |0, 0, — , Oi?) subspace in this case is given by: 

r>12 1 

^cub ^ / , ,, A 12 \ ^;12 , rc C:12 , re ';12 /tin 

-J^ = l^[L^r-i^Z- i^Z,cuh) (^Z + ''+,cubC^+ + <^-XMb<^- (31) 



where 



A 12 

^Z.cub 



m*ujz ( ^x (Cr(-g))' ^ iOY {Cf\-0))' 




— LOy 

2 



wysin (6') L0xcos'^{9) 



m*LOz / exp(±i(/.^)fae^|,Crb(-g) eM±WY)5yci^Uf\-(^) \ .... 
^'^^ " 2^2 ^ ul-Lol ) 

with CT^{-9) = acos(e') + [/? - 7TO*wy/(2?i)] sin(6') and (^"^{-9) ^ -{asm{9) + [f3 - jm*ujx /C^h)] cos{9)). Since 7 
and P have the same sign, the cubic Dresselhaus interaction lessens the contribution of the linear Dresselhaus interaction to 
the effective Rabi frequency under parametric modulation of the electrostatic potential. Similar results are also obtained in 
the case of modulating ujx and tjy; in this case, the effective Rabi frequency is again proportional to the spin-orbit coupling 
parameters only in the presence of a nonzero static electric field, Fx(y) 7^ [otherwise, it is second-order in the spin-orbit 
coupling parameters]. 

Under parametric modulation of both the quan tum d ot's center and oscillator frequency, Rabi oscillations, which depend only 
upon 7, can occur; measuring the Rabi frequency would therefore provide an independent measurement of 7. Consider the 
case when the quantum dot's center is modulated about the y direction at a frequency of ujri, yc{t) — SycSm{ojrit + 4>y), 
while at the same time the quantum dot is being periodically squeezed about the x direction at a frequency of ujr2, ^x{t) = 
ujx + Sllix sm{LUr2t + 4>x)- In this case, efficient spin rotations can be performed in the relevant bimodal Floquet subspace, 
|1) = |0, 0, +, Ipi, 1_F2) and |2) = |0, 0, — , Oj^i, 0^2) when cj^i + Wr2 ~ '^z [here Fx and F2 denote the Floquet states relative 
to the oscillator frequencies w^i and L0r2 respectively]. The effective Hamiltonian in this subspace can be written as: 

f>12 1 

-"cub /',,!,, ,, A 12 a12 a ;fl2 I j12 ;;12 , rl2 -:12 (n-,^ 

-J^ = 2 I'^'-l + ^^^2 -UJZ- i^Z,CMh ~ Oz.cMh) (^Z + t'+,cubC^+ + 0-,cvib<^- (33) 



11 



(A) 



O 



H- 
O 
O 



o 
o 




0.0 



(B) 1.0 



0.68 

time (|j.s) 



1.36 




0.69 

time (fJ-s) 



1.38 



FIG. 4: Numerical simulation of the transition amplitudes, j{0, 0, + |rexp(— i dt'H{t'))\0, 0, — ) p, caused by the simultaneous modulation 
of the quantum dot's oscillator frequency [cjxit) = lux + Slox sin(cij,.2t)] and center [j/c(t) = Si/c sin{LOrit)], (A) with and (B) without the 
Rashba and linear Dresselhaus spin-orbit interaction. Efficient spin transitions, which depend upon the cubic Dresselhaus coupling to first- 
order, can be performed when u>ri + <^r2 ~ ljz [Eq. l[35j]. In both simulations, the following parameters were used: htoy = 0.45 meV, 
hujx = 3 meV, , hujz = 0.1 meV, 6 = 0, SuJx = Wri = 2w^2, Sy^ = 56 nm (equivalent to Ey = 10* eV/m), 7 = 27 x 10"^" eV-m^, 
and (j}x = (t>Y = 0. In (A), a = 4 x 10~^^ eV-m and = 30 nm (giving /? = 2.96 x 10~^^ eV-m), which gave a sizable Bloch-Siegert shift 
of '5z^cub/(27r) = 3.32 MHz. In Fig.|4];A), the observed Rabi frequency was 736.6 kHz, which is within 2% of the calculated theoretical value 



given by Eq. 1 35 1, |i5^^„i,l/(27r) = 721.5 kHz. In (B), the linear Dresselhaus and Rashba spin-orbit coupling constants were artificially set to 
zero, resulting in &^z,cMhl (27i") ~ 400 kHz; the observed Rabi frequency was 724.5 kHz, which is similar to that observed in Fig. 41 A) and is 
closer to the calculated value given by Eq. ( |35^ . This demonstrates that the observed Rabi oscillations are mainly due to the cubic Kesselhaus 
coupling [differences between (A) and (B) are mainly due to higher-order coupling terms involving a and /3]. 



where 

^F.. ^ r^v (34) 

.12 i7cos(6') 5u}x5ycUJri{'m*u}x'^Yf 

The effect of the linear Dresselhaus and Rashba spin-orbit coupling arises mainly in the Bloch-Siegert shift, S'^zcah^ ™ ^Fcub' 
higher-order terms contributions to 5}^^y^^, arising from both the linear Dresselhaus and Rashba spin-orbit coupling can be made 
negligible by using smaller Sujx and 5yc- 

Figure 4 shows the numerical simulation of | (0, 0, +|rexp( — |^ H{t')At')\Q, 0, — ) p under both periodic modulation of the 
quantum dot's center about the y direction and modulation of the oscillator frequency about the x direction in (A) the presence of 
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and in (B) the absence of the Rashba and linear Dresselhaus spin-orbit interactions. In Fig. |4|A), the following parameters were 
used: hujx = 3 meV, hajy = 0.45 meV, hujz = 0.1 meV, dujx = wx/10, Wri = 2wr2, 7 = 27 x 10"^" eV-m^, a = 4 x 10"^^ 
eV-m, — 30 nm (giving (3 « 2.96 x 10"^'^ eV-m), and 5yc = 56 nm (giving Ey ~ 10^ eV/m). With the above parameters, 
'^Fcub/ (2^) ~ 3.33 MHz. The observed Rabi frequency in Fig. Kl A) was 736.6 kHz, which is within 2% of the calculated value 
given by Eq. (35l, |<^±^cubl/(2'''') = 721.5 kHz. Since S^^^^ >1'^±^cubl' ^ri and ujr2 must be tuned to include S^^^^ in order 
to maximize the ampHtude of the Rabi oscillations [in Fig. |4|A), {lUj-i + ujr2 — ^z)/{2tt) = —4.504 MHz+S^^^^J/{2^T)]. Note 
also that in the numerical simulation of Fig. |4|A), the transition amplitude does not go exactly to 1 (maximum value of 0.9861), 
which is most likely due to higher-order corrections to the Bloch-Siegert shift, (Jp^ub- 

Better agreement between the Rabi frequency calculated using Eq. (35i to numerical simulation is obtained in Fig. |4|B) 
where both a and [3 were artificially set to zero, removing the dominant higher-order corrections to both 5^^^^^ and S}^^^^^. In 
Fig. p|B), the observed Rabi frequency was 724.5 kHz, which is closer to the calculated value given by Eq. ( 35 1, | S}^^^^ I / (27r) = 
721.5 kHz. It should be noted that both simulations shown in Fig. |4]give roughly the same Rabi frequency, which, to a good 
approximation, is given by Eq. ( [35| . Thus measuring the Rabi frequency under such parametric modulations should enable 7 to 
be directly determined. It should be noted that for a quantum dot with a non-harmonic potential and in the presence of an out of 
plane magnetic field, it has previousljDIIbeen found that the EDSR technique also generates spin rotations which, to first order, 
depend upon 7 and the cyclotron frequency. 



V. PARAMETRIC MODULATIONS IN SQUARE QUANTUM DOTS 

Besides harmonic potentials, another model electrostatic potential for electrons in lateral quantum dots is that of a square-box 
(hard wall) potential defined by V{x, y) — for < a; < Lx and < y < Ly and V{x, y) — 00 everywhere else. Although 
the hard wall potential is not very realistic in two-dimensional systems, previous studies have utilized such models in analyzing 
transport in quantum dot^SEl and such potentials are the quintessential model for studies of chaos in two-dimensional systems, 
such as in the stadium billiarcPSl. For the square-box potential, the eigenstates for Hq are simply the two-dimensional particle in 
abox states, {r\n,m,±) — ^ -^yj~^sm (^77^^ sin (^^j^^ |±)- In order to induce transitions between the various |n,m,±), 

small, periodic modulations in both the length and the width of the box, Lx and Ly, can be made. Using Eq. Q, the effective 
time-dependent Hamiltonian during modulation of the box's length and width, Lx Lx{t) and Ly Ly{t), can be written 
as: 



2m*L\{t) 2m*L\{t) 2 

^i^^P^ (Ci(-^?)fTx + c,mdz) ~ ^^fy {Ci{e)<Jz + C2i~e)ax) 

dLxit) 1 ^ , A , dLy{t) 1 



dt 2nLxit) (^-^ + ^^-) + ^2^ (^-^ + ^''^ 



where 



Hx = n^|rj, m, ±)(n, m, ±1 

Hy = TO^|n,TO,±)(n,TO,±| (37) 



n.m.ib 



The form of H{t) arises naturally from Eq. previous work on parametric deformations of hard wall potentials have also 
arrived at a similar form for the effective HamiltoniarPEH Note, however, that it is necessary to separate the transformation 
W{t) from the dynamics and only include its effect at the end of the calculatiorP?! 

By performing small modulations of the length and width, Lx Lx + SLx{t) and Ly Ly + 5Ly{t), spin transitions 
can be performed due to the spin-orbit interaction. The effective Hamiltonian in the |1) = |1, 1, +, 1 p) and |2) = |1, 1, — , Op) 
Floquet subspace is given to first-order in dLx{y){t) and for Ur ~ i^z as: 



tt12 

-"eff 



ujz 



12 



(38) 



13 




time (ns) 



FIG. 5: Numerical simulation of the transition amplitudes, |(1, 1, +jTexp(— i J^* dt' H(t'))\l, 1, in a square quantum dot caused by 
either (A) modulating the size of the quantum dot or (B) by applying a time-dependent electric field. In both cases, the Rabi frequency is first- 
order in the spin-orbit coupling strength. The parameters used in both simulations were Lx ~ Ly = 70 nm [giving Ti^tv'^ / (2m*L|-(y)) = 
1.1 meV], Q = 4 X 10"" eV-m, /3 = eV-m, hcuz = 0.1 meV, 9 = 0, (iOr - ujz)/{2tt) = -958.4 kHz, and (j}x = 0. In (A), 
Lxit) = Lx + SLx smiujrt) with 5Lx = 5 nm, which gave an effective Rabi frequency of 10.6 MHz, which is in excellent agreement with 
the calculated value given by Eq. i40i, \5±\/{2'it) = 10.6 MHz. In (B), the apphed electric field was E{t) = Ex s\n{ijjri)x with Ex = 10** 
eV/m, which gave an effective Rabi frequency of 3.94 MHz, which is in excellent agreement with the calculated value given by Eq. (|42j, 
\5f\/{2Ti) = 3.92 MHz. 



where 



^ 32m* ^ n^l + (-!)") / (CiM)f^x , {U^O)f m \ 

^ 4y n^(l + (-!)") / Ci{-0)vxSLxeM±i<f>x) _ C2{-0)vySLy eM^^M 

(40) 

where rjx(Y) = '2,m*L0zL'^x(Y)li'^''^'^)- choosing LOr = ioz — ^^z^ efficient spin transitions can be generated in square quan- 
tum dots [Fig. [sJA)] since the Rabi frequency, |(5± |/(27r), is directly proportional to the spin-orbit coupling strength, in contrast 
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to a parabolic quantum dot undergoing modulations of confining frequency and in the absence of a static electric field [Eq. ( 23 1] . 
Figure|5|A) presents an exact numerical simulation of the transition amplitude, 1, +|rexp(— i J* dt' H{t') /h)\l, 1, — )p, in 
the presence of modulating the effective size of a square quantum dot. The following parameters were used in the simulation: 
Lx ^ Ly ^ 70 nm [giving H^t:^ /{2m* Lj^^y)) ~ 1-1 and tjx{y) = 0.087], SLx = 5 nm, a = 4 x lO^^^ gV-m, (3^0 
eV-m, huz — 0.1 meV, 9 — 0, and (w,. — ujz) /(27r) = —958.4 kHz. Such parameters gave an effective Rabi frequency of 10.6 



MHz, which is in excellent agreement with the calculated value given by Eq. (40 1, |(5+ |/(27r) — 10.6 MHz. 

Alternatively, an electric field, E{t) — Ex sm{ujrt + (j>x)x + Ey sm{ujrt + (j)Y)y, applied to the quantum dot can also induce 
single spin rotations by using EDSR effects. Incorporating the interaction with the electric field, V{t) — —eE{t) ■ r, an effective 
Hamiltonian in the |1) and |2) subspace [when iVr ~ ojz] can be written as: 



tt12 

^EFF _ 1 

n ~ 2 

where 

,2 



Ur -UJZ+ df + Jfal^ + ^EF^l^Z (41) 



^EF ^ 32m;^ + / Lj,ExCi{-9)iix exp(±i0x) _ Lj- Ey C2{- 9)riY (^M±^<Py) 



(42) 



Like the case of modulating the size of the quantum dot, the effective Rabi frequency, |(5^E|/(27r), is again propor- 
tional to the spin-orbit coupling strength. Figure |5jB) presents an exact numerical simulation of the transition amplitude, 
1(1, 1, +|rexp(-i /J dt'H{t')/h)\l, 1, in the presence of an electric field, E{t) = Ex sin(cjri + (j)x)x. The follow- 
ing parameters were used in the simulation: Lx = Ly — 70 nm [giving h' t:^ / (2m* L"^^ f^y)) ~ 1-1 meV and ?7x(y) — 0.087], 
Ex = 10^ eV/m, a = 4 x lO^^^ eV-m, /3 = eV-m, huJz = 0.1 meV, 9^0, and {ur - t^z)/(27r) = -958.4 kHz. Such 
parameters gave an effective Rabi frequency of 3.94 MHz, which is in excellent agreement with the calculated value given by 



Eq. (42 1, \df\/{2Tr) = 3.92 MHz. 



VI. PARAMETRIC ORBITAL AND SPIN EXCITATIONS FOR INCREASING SPIN POLARIZATION 

In addition to single spin manipulations in parabolic quantum dots, combined transitions between both the spin and the orbital 
degrees of freedom can be generated under parametric modulation of V{X, Y). In a parabolic quantum dot, application of an 



electric field [H{t) in Eq. (25 i] and/or modulation of ux and ujy [H{t) in Eq. (20l] can directly induce both orbital and spin 
excitations. Such transitions are an essential component for proposals to increase spin polarization in a quantum dot's ground 
electronic state'*". Consider first the case of parametrically modulating a parabolic quantum dot's oscillator frequency, ujx, in 



the absence of any static electric fields [H{t) in in Eq. (20i with Fx = Fy = eV/m]. In the following,we will neglect the 



cubic Dresselhaus interaction [i/™'' in Eq. (|30|l] and will take uy > uix in order that the parametric modulation can be chosen 



to connect states involving changes in only one of the orbital degrees of freedom, such as |0, 0) going to |1, 0). If we wish to 
generate transitions between the states 1 1, 0, +) and |0, 0, — ), ujx must be modulated at frequency of ujr ~ llix — ^z- Labeling 
the relevant states as |2) |0, 0, — , If) and |3) = |1, 0, +, Of), the effective Hamiltonian in the {|2), |3)} Floquet subspace 
under parametric modulation of ux (at uJr « ujx — ujz) is given by: 

&23 1 

-"EFF ,, , ,, , A 23 , r23\ -;23 , r23'-23 , c23':23 ,41^ 

= 2 i - UJX + + i-^Z + "z ) ^Z + "+^+ + ^- (43) 



A 



23 



m*LOz (2uJx{U-0)f , u;Y{C2{-0)f 



S 



23 _ f i^X - ^zf 

^ Ildx \4w|- - {ijJx - (^zY 



4 = \h;r-^ 2 2 exp ±i(/.x (44) 

By choosing ojr — lox — loz — ^z ^ ^^z^ efficient spin and orbital transitions between the states |2) and |3) can be generated 
[Fig. pi A)] since the Rabi frequency, |(5|f |, is directly proportional to the spin-orbit coupling. Note that there can be a significant 
Blocl>Siegert shift, (5|?, since the transition involves different electronic orbitals. 
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Figure|6|A) presents an exact numerical simulation of the transition amplitudes, 0, +|Texp(— i J* dt'H{t') /?i)|0, 0, — )p 
(black curve) and |(0, 0, +|Texp(— i dt' H{t')/'h)\0,0,+)\'^ (red curve), in the presence of the modulation u>x{t) — uJx + 
Sujx sm{LUrt). The following parameters were used in the simulation: Tiujy = 1 meV, fiojx = 0.25 meV, a = 4 x 10^^'^ eV-m, 
P = eV-m, hcvz = 0-1 meV, 9 — 0, Fx — Fy — eV/m, and 5ujx — "^x/lO- Such parameters gave an effective Rabi 



frequency 48.1 MHz, which is in excellent agreement with the calculated value given by Eq. (44 1, |(5^'^|/(27r) = 48.3 MHz, 
under "on-resonant" modulation [(w^ - (wx - t^z))/(27r) = 32.32 MHz -5f /{2tt) = 17.35 MHzf Note that the state |0, 0, +) 
essentially does not evolve into any of the other states [the red line in Fig. |6|A)] thus confirming that the electron must be spin 
down in order to be excited from the state |0, 0). As mentioned earlier, such spin selective transitions can be used to increase the 
spin polarization of a quantum dol"**^, which will be discussed later in this section. 

Since the relevant transitions involve both orbital and spin excitation, the required modulation frequency is quite large, ujr ~ 
ujx — ~ 36.3 GHz, which can be experimentally challenging. However, higher-order processes which utilize smaller 
values of lUj. can be found in order to induce efficient spin and orbital transitions. The simplest higher-order process to induce 
transitions between the states |1, 0, +) and |0, 0, — ) occurs when ~ {^x — '^z)/2. In this case, the Floquet states |1, 0, + , 0^) 
and |0, 0, — , 2i?) are degenerate. The effective Hamiltonian in this subspace is given by [for lOj. k, (ojx — loz)/'^]'- 

f>23 1 

2ti; ^ / , ,, , o, , , A 23 , c23 N ;;23 , r23 ^^23 , c23 ^^23 (ac\ 



where 



2 



5uj\ ( {'^x - ^zY 



c23 ^ "'^X 



^2t,± = ^^h;, ^7 ^1 T ho ^exp(±i20x) (46) 

In this case, the effective Rabi frequency is again directly proportional to the spin-orbit coupling strength and is also proportional 
to the square of the modulation strength, {SuixY- 

Figure|6|B) presents an exact numerical simulation of the transition amplitudes, 0, +|Texp(— i At' H{t') /h)\Q, 0, — )p 

(black curve) and | (0, 0, +|rexp(-i g&t' H{t') in)\Q,Q,+)\'^ (red curve), in the presence of modulating the oscillator fre- 
quency of the parabolic well at ujr ~ {(^x — '^z)/2. The following parameters were used in the simulation: fiLOy — 1 
meV, n^x = 0.25 meV, a = 4 x IQ-^^ gV-m, /3 = eV-m, Uujz = 0.1 meV, 6* = 0, Jwx = wx/10, and 
{ujr - l/2(t^x - ujz))l^'^) = -Sjl z/(47r) + 16.16 MHz = 14.42 MHz. The above parameters gave an effective Rabi 
frequency of 8.77 MHz, which is within 6% of the calculated value given by Eq. (46 1, |(5i^|/(27r) w 8.26 MHz. Note also that 



the state |0, 0, +) doesn't evolve during the parametric modulation of ujx [red curve m Fig. |6|B)]. 

In addition to modulating the confinement frequency, the spin-orbit coupling constants, a and /3, can also be modulated in 
order to produce combined spin and orbital excitations which are first-order in a and [3. By modulating an electric field generated 
by a surface gate above the quantum dot, the Rashba spin-orbit coupling parameter within the quantum dot can be controlled '3 
additionally, for quantum dots formed using semiconductor heterostructures, surface gates could also be used to change 
thereby changing the linear Dresselhaus coupling constant, /3 = 7(P|). Such modulations of the spin-orbit interaction have 
been previously suggested as a way to perform both spin and orbital excitations in the absence of a magnetic field'"*. Finally, 
it should be noted that application of an electric field can only couple the states |1, 0, +) and |0, 0, — ) to second-order in the 
spin-orbit coupling for a parabolic quantum dot; therefore, the EDSR technique cannot generate such transitions as efficiently as 
modulating ojx when the spin-orbit coupling is weak. 



A. Increasing spin polarization in tlie lowest electronic oribtal 

Throughout this work, relaxation has been neglected when calculating the various spin transitions generated under parametric 
modulation of the quantum dot's electrostatic potential. In the previous sections where only effective spin rotations within the 
lowest electron state were considered, neglect of both Ti and T2 relaxation se emed just ified since the calculated Rabi frequencies 
(1 — 30 MHz) were one to two orders of magnitude larger than the measurecP^Eulllol x /y^ values in quantum dots (we assume 
that I/T2 is on the same order"" as l/Ti). However, relaxation of excited electronic states in quantum dots occurs on the 
nanosecond time scal^, which is much faster than the calculated Rabi frequencies (order of 10 MHz) associated with combined 
spin and electronic excitation shown in Fig. [6} therefore, relaxation effects must be considered when examining such transitions. 

The ability to coherently couple the states |0, 0, — ) and 1 1, 0, +), coupled with relaxation, can be used to help spin polarize the 
electronic ground state of a quantum dot as depicted in Fig. |7jA). From Fig. |6] parametric modulation of the oscillator frequency 
can connect the state |0, 0, — ) to the state |1, 0, +) while leaving the state |0, 0, +) relatively uncoupled from any other state of 
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FIG. 6: (Color online) Numerical simulation of the transition amplitudes, |(0, 0, — |Texp(— i dt' H(t'))\l,0, +)f (black curve) and 
I (0, 0, +|r exp(— i df'_ff (t'))|0, 0, +)p (red curve) caused by modulating the oscillator strength defining the parabolic quantum dot, when 
(A) LJr ~ LOx — ijJz and (B) ujr ~ (<^x ^ ^z)/'2,. The following parameters where used in the simulations: Tiujy = 1 meV, hcJx = 0.25 
meV, a = 4 X 10~^^ eV-m, P — eV-m, huz = 0.1 meV, 6 — 0, Sujx = and 0x = 0. In (A), modulation of ujx at a frequency 
of ujr/(2n) ~ (ujx — ioz)/{2n) = 36.3 GHz can generate effective transitions between the states |0, 0, — ) and |1,0,+), which gave a 
Rabi frequency of 48.1 MHz, which is in excellent agreement with the calculated value given by Eq. \5f\/{2Tv) = 48.3 MHz. In (B), 
modulation of lux at a frequency of ujr/(2n) ~ {ujx — a;z)/(47r) = 18.15 GHz resulted in a Rabi frequency of 8.77 MHz, which is within 
6% of the calculated value given in Eq. l |46[ l, |52E.±l/(27i") ~ 8.26 MHz. Although the resulting Rabi frequency is smaller in Fig. [6|B) than 
in Fig. |6jA), the required ujr is a factor of two smaller, which makes such modulations more experimentally feasible. 

the quantum dot. An electron in the excited state |1, 0, +) can quickly relax to the state |0, 0, +) due to the direct coupling of 
the electron to piezo-phonons^^, whereas the relaxation of the state |1, 0, +) to the state |0, 0, — ), which requires an effective 
spin-phonon coupling, occurs at a much smaller rate. This differential relaxation can be used to generate spin polarization in the 
lowest electronic state in a manner similar to optical pumping"*"* and Overhausei'*"* techniques. 

In the following work, the dynamics under parametric modulation of the confining potential has been restricted to the 
following subspace for simplicity: {|1) = |0,0, +), |2) = |0, 0,— ), |3) = |1,0,+)}, and the calculations were per- 
formed using different values of S'^ taken from the Fig. |6] The density matrix in this subspace can be written as p(t) — 

+ P23(i)|2)(3| + P32(i)|3)(2| (where only coherence between the states |2) and |3) have been considered). 
Defining Wij to be the transition rate from state i to state j, and to be the decoherence rate for the coherence between states 
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1 2) and |3), the various coefficients in p{t) can be found by solving: 
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(47) 



where the parametric modulation of the quantum dot is assumed not to affect the values of Wij and r23 in Eq. (47 1. In the 
absence of any term connecting the states |2) and |3), i.e., 5'^ — 0, the various transition rates, Wij, must satisfy the following 
conditions in order to ensure that the equilibrium density matrix, peq = Si Pii I*) (*|, is a solution to Eq. (47 1: 



eq ''^ 



where is the equilibrium population for state j, with p^^^ > P22 > P33 J^iPu — 1- 

In the presence of nonzero 6"^, the effective spin polarization of the ground electronic state, Pz 



,eq 



its equilibrium value, 



P11+P22 



(48) 



, can increase from 



-P22 



1-P22 



Such an increase occurs since W31 involves only orbital relaxation, which is in general 



much faster than the transition rate 14^32, which requires both orbital relaxation and a spin flip. Applying parametric modulations 
to the quantum dot can therefore transfer population from state |2) to state |3), which subsequently relaxes to state |1), thereby 
increasing the population difference between states |1) and |2). This is shown in Fig. [tJA). 

The steady-state spin polarization under such parametric oscillations can be found by setting the left-hand side of Eq. ( [47] i to 
zero and solving for pn and P22, which gives: 



P 



Steady-State 



W31 (Wi3 + W23 + + (1^21 - W^2 - W^13) (m^31 + W^32 + 2^^) 

W31 (m^23 - VF13 + 2^) + (1^21 + Wx2 + W13) {W^x + W32 + 2^) 



(49) 



Alls 



The transition rate, W31, is mainly dominated by the coupling of the electron to piezo-phonons. From Eq. (4) of 
rate is given as Wzi ~ 3.74 x 10* s^^ (for htox — 0.25 meV and Tiuy — 1 meV); using Eq. (7) of Ref."^^ as a rough estimate 
for W32, one obtains W32 ~ 5 x 10^2 g-i^ which satisfies W31 ^ M^32- Although no theory has been attempted in this work 
to calculate r23, as long as 2r23 > VF31, the diagonal elements, pu, will be nonnegative (as seen from numerically integrating 



Eq. (|47|i). Although the value of r23 does not drastically affect p^"='"'5''''''"'=^ \[ (joes affect the time scale in which 
is reacEed. In the following we simply take, for illustrative pu rposes, r23 = IOIV31. Finally, since the effective Ti for GaAs 
quantum dots is on the order of milliseconds to microseconds '^^UlSlZOI ^ conservative value of Ti = 100 /is was chosen, which 
provided the values of W12 and W21 from the condition that l/Ti = W12 + W21 — W2i(l + cx-p{TiuJz /{ksT))) [where the 
populations were assumed to be given by the Boltzmann distribution^. For T = 2 K and hojz = 0.1 meV, W21 ~ 3.59 x 10^ 
s~^. Using the parameters from Fig. pfA) (|(5^3|/(27r) = 48.1 MHz), Pz{t) was found by integrating Eq. (47 1 and is shown 
(solid curve) in Fig. RB), with p^tea(j-state q g2 ~ 2.82P^'', which is reached within 0.15 /is. No oscillations are present in 
Pzit) since T23->\df\. 

For the case of modulating utx 31 a frequency of ujr ~ (wjc — uJz)/'2' [with |(5|3 _|_ |/(27r) « 8.77 MHz as shown in Fig. pfB)], 



it took roughly twenty times as long (w 3/is) to reach P< 
It should be noted that p^"=''''5'"'''^"= 
-P33 



steady-state 



0.62, which is shown by the dashed curve in Fig. 7[B). 
is approximately equal to the initial equilibrium population difference between the states |1) 

cq 

and 13), /"gp « 0.62. Thus turning on the coupling between 12) and 13) allows the initial "orbital" polarization between the 

states |1) and |3) to be transferred into spin polarization between the states |1) and |2). This process is similar to the Overhauser 
effect, where, by coupling the electron and nuclear spins, the initial electron spin polarization can be converted into nuclear 
polarization'*''. Therefore, in order to increase p^"='"*>'''*"'"=^ experiments should be performed at low temperatures, large magnetic 
fields, and with increased lateral confinement, i.e., large ujx- However, increasing ujx requires modulating the quantum dot at 
a higher ojr (which may be unfeasible experimentally) in addition to the fact that \S'^\ decreases with increasing ojx, thereby 
increasing the time it takes to reach p^''=^<iy-''''"'= 

Finally, it should be noted that the above calculations assumed that the quantum dot was isolated and absolutely closed from 
the leads, i.e., no electrons could enter or exit from the dot. However, if electrons are able to tunnel in/out of the quantum dot, 
modulating the electrostatic potential could be used to selectively "kick" out spin down electrons. "Re-zeroing" the state of the 
quantum dot to spin up could be used as an initialization step for a possible quantum computatior(i2li21 Such a process would 
require the tunneling rate out of the dot for an electron in state |3) to be much faster than W31. However, the evolution under 




FIG. 7: Method for increasing the spin polarization of the lowest electronic state, |0, 0). In (A), parametric modulations are performed in 
order to coherently connect the states |2) = |0, 0, — ) to the state j3) = |1, 0, +) (which is denoted by Sji^), thereby transferring population 
from state |2) to state |3). The population in state |3) predominately relaxes to the state jl) = |0, 0, +) due to the electron directly coupling 
to piezo-phonons, thereby increasing the spin polarization of the ground electronic state, Pz ~ p\\+p22 • Integration of Eq. |47| gave 
Pzit) for parametric modulations of ujx at the frequency ujr ~ — (solid curve) and at l^It = (i^x — uJz^j'i- (dashed curve). Along 

" 1^31 = 3.74 X 10« s-\ W32 = 5 X 10-2 

8.77 MHz (dashed curve). 



iteady-state 



2 K, Taa = lOW^si, and i<5±'l/(27r) 
0.62 [Eq. ([49}], although for larger \&\ 



with the parameters used in Fig. [6] the following parameters were used in solving Eq. 1 47 
W21 = 3.59 X 10^ s-\r 
In both cases, P\ 



48.1 MHz (solid curve) and |5|^|/(27r) 



psteady-state 



was reached on a faster time scale [0.1 5^s (solid curve) 



compared to 3pis (dashed curve)]. 



parametric modulations of the confining potential would have to be reconsidered, since in the case of a non-isolated quantum 
dot, there would be a finite probability during the parametric modulation that an electron in the state |1) would get excited and 
kicked out of the dot if the tunneling rate out of the dot is very large. 



VII. DISCUSSION AND CONCLUSIONS 



In this work, a general formalism combining Floquet theory with effective Hamiltonian theory was used to study the spin dy- 
namics under parametric modulations of a lateral quantum dot's electrostatic potential in the presence of spin-orbit coupling. In 
parabolic quantum dots, modulating the center of the quantum dot, i.e., performing an EDSR experiment, was found to generate 
larger Rabi frequencies than parametrically modulating the confining frequency (in the absence of static electric fields), since 
the latter is second-order in spin-orbit coupling. However, in the presence of a static electric field, both methods gave similar 
Rabi frequencies for a parabolic quantum dot. For square dots, both EDSR and modulating the width/length of the quantum dot 
generate Rabi frequencies which are first-order in spin-orbit coupling. The modulation frequency must be chosen with precision 
on the order of the Rabi frequency (on the order of tens of MHz) in order to maximize the amplitude of the Rabi oscillations, 
thereby providing better spectroscopic precision of the quantum dot's energy levels relative to transport measurements, where 
thermal effects decrease the spectral resolution. 

Inclusion of the cubic Dresselhaus spin-orbit coupling didn't dramatically alter the results obtained for parametric modula- 
tion in a parabolic quantum dot. However, the cubic Dresselhaus interaction tends to decrease the contribution of the linear 
Dresselhaus interaction to the Rabi frequency. Furthermore, measurement of the Rabi frequency under, for example, an EDSR 
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experiment using different orientations of the in-plane magnetic field could be used to determine the ratio of the Rashba spin- 
orbit coupling constant to the difference of the linear Dresselhaus spin-orbit coupling constant with the product of the cubic 



Dresselhaus coupling constant and the oscillator frequency {ojx or ojy). For example, using Eq. (29 1 and including the cubic 



Dresselhaus coupling from Eq. (32 1 gives: 



ZZ^-) . ^—^ (50) 

In order to separate 7 from (3 in Eq. ( [50| , experiments could be repeated for different values of cuy, which would leave the linear 
Dresselhaus's contribution to ZZ{Tr/2) unchanged but would alter the cubic Dresselhaus's contribution. Alternatively, it was 
shown that applying a time-dependent electric field along the y direction coupled with parametric modulation of the confinement 
frequency along the x direction could be used to generate Rabi frequencies which were proportional to 7 and were independent 
of a and to first-order Measurement of the resulting Rabi frequency would provide another independent measure of the cubic 
Dresselhaus spin-orbit coupling constant. 

Parametric modulations of the confining potential were also shown to generate combined spin and orbital excitations in a 
parabolic quantum dot with Rabi frequencies on the order of tens of MHz. However, since electronic relaxation times are on the 
order of nanoseconds, the effects of orbital relaxation had to be taken into account. A combination of coherent spin and orbital 
excitation with orbital relaxation was shown to be able to increase the spin polarization of the ground electronic state [Fig. |7| 
by transferring the initial "electronic polarization" between the ground and excited electronic state into spin polarization of the 
ground state. For the parameters chosen in this work, this corresponded to a three-fold increase in the spin polarization of the 
ground electronic state. Larger spin polarizations could be achieved by either using a more confined quantum dot (larger cox), 
larger magnetic fields, or by going to lower temperatures. 

Finally, the spin control developed in this work used only parametric modulations of idealized electrostatic potentials 
(parabolic and square-box), which should only be considered as a model for spin control in lateral quantum dots. For such 
methods to be used in actual experimental quantum dots, more realistic electrostatic potential^ for a quantum dot, i.e., non- 
parabolic and non-square-box potentials, are required. Furthermore, additional work is needed in order to better characterize the 
actual time-dependent electrostatic potentials generated by modulating the surface gate voltages of the quantum dot along with 
designing optimal configurations of surface gates in order to generate a desired transition. The results presented in this work 
could furthermore be extended to many-electron quantum dots, where the effects of electron-electron coupling on perform- 
ing spin excitations should be examined. Recent theoretical worlP^ljias demonstrated that electronically controlled, magnetic 
dipole-like couplings between spins in different quantum dots can also be generated. These couplings, along with single spin 
rotations, could be used to generate multiple-quantum coherences between the electrons in different quantum dots by completely 
electrostatic means. In addition to being used to better characterize the electron-environment coupling of such highly correlated 
states, generating such highly correlated states could be used as an intermediary step during a quantum computation. 
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APPENDIX A: THE EFFECTIVE HAMILTONIAN 

If for a given Hamiltonian (written for convenience in Floquet space), Hp — Hp + Vp, there exists a subspace of interest, 
Q, which is weakly coupled by Vp to the rest of the Floquet space, U, then the dynamics within Q can be separated from U by 
constructing an effective Hamiltonian from Hp, H^^, which is block-diagonal in the Q and U subspaces^^ Defining Pq to 
be the projection operator onto the subspace Q, Pq ~ J2\a m.F)eQ 1"^' ^f) i^^ '^t^fI, and P[/ to be the complementary projection 
operator onto the subspace U, Pjj ^ Ip — Pq, H^ can be determined by constructing a unitary transformation, exp{Sp), such 
that exp (^Sp^ Hp exp (^—Sp^ = H^^, where in order to ensure that iJfP^ is block diagonal in both the Q and U subspaces, 
H^^ must satisfy the following: 

pQHf^Pu = PuHf^Pq - (Al) 

Using Eq. ( Al i, a perturbation expansion for Sp in powers of Vp, Sp = J2m=i '^f"'*' constructed. Separating Vp into 

its diagonal and off-diagonal components, Vp = Vp+Vf, with Vp = PqVp Pu + PuVp Pq and Vf = PuVpPu + PqVpPq, 
Sp can be written up to m = 3 as: 



/ ic(i)ifl \ {a,mp\V^\P,np 

{a,mp\S'p'\P,np) - 



Ea ~ + 'h{m ~ n)ujr 
(2),. , _ {a,mp\[s'p'\vf]\(3,np) 



{a,mp\S},'\P,np, - ^ , n 

ha — hp + rL[m — n)ujr 

I 1^(3),. . _ {a,mp\[sf ,VP]%np) {a,mp\[S^P ,[S^P ,VmP,np) 

'^a,mp\Sip \p,np) — — - r \ — — = ^, ^ ^ (a/; 

Ea — Ep + n[m — n)ujr A {ha — hp + nym — n)ujr) 
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Using the above values of Sp in Eq. ( A2 1, the effective Hamiltonian in the subspace Q, PqH^Pq — J2 
(up to order (Vp)^) as: 
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